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Abstract 



'■^ ■ Properties of the two-dimensional ring and three-dimensional infinitely 

C3 ' long straight hollow waveguide with unit width and inner radius po in the 

O , superposition of the longitudinal uniform magnetic field B and Aharonov- 

, >^/ Bohm flux are analyzed within the framework of the scalar Helmholtz 

equation under the assumption that the Robin boundary conditions at 

the inner and outer confining walls contain extrapolation lengths Ai„ and 

^ ' Aout, respectively, with nonzero imaginary parts. It is shown that, com- 

^\ I pared to the disk geometry, the annulus opens up additional possibilities 

00 . of varying magnetization and currents by tuning imaginary components of 

C^ ' the Robin parameters on each confining circumference; in particular, the 

possibility of restoring a lossless longitudinal flux by zeroing imaginary 
[■*-. ■ part Ei of the total transverse energy E is discussed. The energy E turns 

(.~:> ' real under special correlation between the imaginary parts of Ai„ and Aout 

with the opposite signs what corresponds to the equal transverse fluxes 
through the inner and outer interfaces of the annulus. In the asymptotic 
case of the very large radius, simple expressions are derived and applied 
to the analysis of the dependence of the real energy E on Ai„ and Aout. 
^\^ ' New features also emerge in the magnetic field influence; for example, if, 

JH , for the quantum disk, the imaginary energy Ei is quenched by the strong 

intensities B, then for the annulus this takes place only when the inner 
Robin distance Ain is real; otherwise, it almost quadratically depends on 
B with the corresponding enhancement of the reactive scattering. Closely 
related problem of the hole in the otherwise uniform medium is also ad- 
dressed for real and complex extrapolation lengths with the emphasis on 
the comparative analysis with its dot counterpart. 



1 Introduction 

Recent analysis [1] extended into the whole complex plane earlier theoretical 
research [2-4] on the influence of the real Robin parameter A [5] on the properties 
of spatially confined two- (2D) and three-dimensional (3D) domain il. Obtained 
dependencies on the example of the infinitely long cylinder and its circular cross- 
section generalized previously known results for different real A without [6-11] 
or with [7, 12-24] axial magnetic field B. Robin parameter also known as the 
extrapolation length appears in the boundary condition at the sample interface 
dn with its inward normal unit vector n 



nV*la,, - j^ 



(1) 



for the function ^(r) being a solution of the scalar Hclmholtz wave equation 

V2*(r) + fc2*(r) = 0. (2) 

Here, the wave vector k is equal to the ratio of the frequency of the electro- 
magnetic or acoustical oscillations uj to the corresponding speed of propagation 

.^"- (3) 

while for the quantum mechanical particle with mass rrip it is expressed as 



k = ^2mpE/h, (4) 

where i<^ is a total energy of the particle and h is the reduced Planck constant. 
In the former case, a solution ^(r) is an acoustical or electromagnetic potential 
through which an air pressure and velocity or intensities of the electric and 
magnetic fields are expressed while in the latter case it is a wave function with 
its square defining the probability of the location of the particle and being 
normalized according to 

/" l*(r)l^dr = l. (5) 

The same boundary condition (1) applies also [25-29] to the solution of the 
nonlinear Ginzburg-Landau (GL) equation [30] 

V2*(r) + fc2*(r)+/31*(r)l^*(r) = (6) 

with positive GL parameter /3. Despite its venerable age, phenomenological GL 
theory of superconductivity [30,31] continues to be a powerful tool of studying 
superconductors as its predictions are in a very good agreement with the exper- 
iment. Order parameter ^(r) from (6) defines the density of superconducting 
particles rig, 

ns - |*(r)|^ (7) 

while the role of the energy E is played by the GL parameter —a which is 
expressed via the actual temperature T of the superconducting material, the 



bulk critical temperature at zero magnetic field Tc, zero-temperature coherence 
length ^ (0) and Cooper pair mass being equal to the double bare electron mass 

TOe: 

Thus, minimizing the lowest eigenenergy leads to the increase of the critical 
temperature T. 

It is "worthwhile to note that de Genncs distance, as the extrapolation length 
A is called in the language of this field of scientific research, can take negative 
values for the border with the other superconductor with higher critical temper- 
ature what leads to the increase of its own T. This theoretically predicted en- 
hancement of superconductivity [11,32-34] was indeed observed in cold worked 
Ino.993Bio.007 foils [32] and tin samples [35,36]. 

If a magnetic field B is applied to the system, its infiuence is included into 
Eqs. (2) and (6) via the vector potential A: 

'-^T-^) *(r) + fc2^(r) = 0, (9) 

jIa) *(r) + fc2*(r) + /3|*(r)|^*(r) =0, (10) 

where q is the particle charge; in particular, for the Cooper pair, q = — 2e, with 
e being the absolute value of the electronic charge. Here, 

B = VxA. (11) 

Vector potential enters also the GL boundary condition: 



h 



dn A 



(12) 

an 

In the GL theory [30, 31], Eq. (10) with boundary condition (12) has to be solved 
together with an expression for the supercurrent density js 

j, =i-^(*V**-**V*)-— A|*|^ (13) 

2mp vrip 

and the Maxwell equation 

V X B = /ij„ (14) 

/x is a magnetic permeability. System of equations (10), (13), (14) was analyzed 
in the Neumann limit, A = 00, of condition (12) for the solid cylinder [37- 
44] and for studied below geometry of its hollow counterpart [45,46] while the 
calculations of the combined infiuence of the magnetic field and real de Genncs 
distance based on either linear, Eq. (9), [16,17,19,24,47-55] or nonlinear, Eq. 
(10), [56-64] GL theory revealed, for different shapes, a significant infiuence of 
the parameter A on the nucleation of the superconductivity, critical magnetic 
fields and localization properties of the superconducting state. We remark that 



the linearised GL equation (9) which is the main subject of the present study 
correctly captures the physical phenomena in the uniform magnetic field B and 
temperature T ranges close to the transition to the normal state when the order 
parameter ^(r) is small and, accordingly, the cubic term in (10) can be safely 
neglected [31]. The influence of the cubic term can be strongly suppressed 
also by the choice of the metal or alloy [65] since the GL parameter j3 contains 
material-dependent density of states at the Fermi energy A^(0), coherence length 
^(0), critical temperature Tc and the mean free path I [31]. In addition, it can 
be shown that equations derived below for the linear case follow also from the 
complete nonlinear GL theory [52,53]. Moreover, a comparison between the 
theory [66] and experiment [67] revealed that, in the Neumann limit, A = oo, 
the linearized GL equation correctly captures features of the aluminum disks 
of the different radii when the predictions of the full GL theory sometimes are 
in a worse agreement with the experiment [67] than the results produced by 
the solution of Eq. (9). Thus, the analysis of the linear magnetic Helmholtz 
equation is indispensable in the study of the properties of superconductors. 

If the extrapolation length A can take negative values, a natural question 
arises: what happens if the Robin parameter is complex! Attempts to answer it 
have been made during the study of the scattering phenomena in different media 
[68] such as a sound duct with porous lining [69-76] , impedance electromagnetic 
waveguides [77] , ferrite- filled resonator systems [78] , absorbers in high-frequency 
electromagnetic scattering [79]. A comprehensive answer showed that for the 
infinitely long cylinder with singly connected circular cross section the imaginary 
part of the transverse complex energy E^ exhibits a pronounced maximum as a 
function of the imaginary part Ai of the de Gennes distance [1] . The energy E^ 
is a transverse component of the total energy E entering equation (2) via (4). 
As a result, the current undergoes a resonant alteration as it flows down the wire 
exponentially increasing or decreasing (depending on the sign of A^) with the 
longitudinal distance. This change of the longitudinal flux is accompanied by 
the transverse radial currents through the circumference of the disk absent for 
the real de Gennes distances [31]. Nonzero real part A^ of the Robin parameter 
and axial uniform magnetic field B quench the resonance and, for their large 
values, restore the lossless longitudinal current. Physically, it is explained by 
the fact that, for example, the increasing magnetic field squeezes the charged 
particle to the cylinder axis; accordingly, with magnetic intensity growing the 
influence of the boundary decreases and transverse energies E-^ transform, for 
the large B, into the Landau levels. 

Apart from mentioned above acoustical and electromagnetic systems, com- 
plex extrapolation lengths can be realized for the superconducting materials 
too [1]. Namely, it was argued recently [80-82] that the electric field £ applied 
perpendicularly to the surface should be accounted for in the total de Gennes 
distance Atot by the addition to the inverse zero- field extrapolation length 1/A 
of the extra term £/Us 

iUot A g^Q Us 



with the potential Us being expressed through the parameters of the GL theory: 

l^^,|ln7;Mi^ (16) 

Here, dimensionless GL parameter n is the ratio of the zero-temperature London 
penetration length A(0) to the coherence length, k = A(0)/^(0); e^ is supercon- 
ductor ionic background permittivity, and c is speed of light. If the permittivity 
e^ has a noticeable imaginary part, so does, according to (15) and (16), the total 
de Gennes distance too. In this model, a variation of the imaginary part of the 
Robin parameter is achieved by a simple change of an applied gate voltage. 

Considered in Ref. [1] model of the solid cylinder offers only one surface 
through which the transverse flux can enter or leave the sample. One can ex- 
pect that the multiply connected structure with at least two confining interfaces 
offers additional channel(s) of controlling its transport and thermodynamical 
properties. The main subject of the present research is to investigate the inter- 
action of the inner and outer complex Robin parameters and its influence on 
the conductivity and magnetization of the hollow cylindrical waveguide and its 
2D annular cross section. It is shown that, indeed, varying the signs and mag- 
nitudes of the de Gennes distances on each confining wall, one can manipulate 
these properties in a wide range; in particular, the conditions for the lossless 
longitudinal flow down the channel are derived and analyzed: if the incoming 
transverse torrent through one surface with complex A is equal in magnitude to 
the outgoing flow through the other wall with the opposite sign of the imaginary 
part of its extrapolation length, then the longitudinal current does not change 
along the duct. New features emerge also when the magnetic field is applied 
parallel to the channel axis; namely, depending on the sign and magnitude of 
the inner de Gennes distance, the alteration of the current can be manipulated 
at will by the appropriate change of B. For example, the lossless current down 
the waveguide at the strong magnetic intensities is asymptotically achieved only 
for the real inner Robin parameter; otherwise, since the field pushes the charged 
carrier closer to the reactively scattering inner surface, the attenuation or am- 
plification of the longitudinal flux almost quadratically depends on the field. As 
a by-product, not considered before aspects of the problem of the hole in the 
film are discussed too and its comparative analysis with the solid cylinder case 
is performed. 

The paper is organized as follows. In Section II our model is presented 
and a necessary formulation of our method is given. Section III is devoted to 
the presentation and detailed mathematical and physical interpretation of the 
calculated results. Summary of the research is provided in Section IV. 

2 Model and formulation 

Infinitely long 3D straight wire of the annular cross section with the inner po 
and outer (po + d) radii is placed into the uniform magnetic field B with its 
direction coinciding with the channel axis (Fig. 1). For the zero inner length. 



Figure 1: Schematic of the infinitely long hollow superconductor wire of the 
inner radius po a-nd the outer one (po + d) subjected to the superposition of 
the uniform magnetic field B and the AB flux $ parallel to its axis. Channel 
walls support Robin boundary condition, Eq. (1), with the complex outer Aout 
and inner Ai„ extrapolation lengths uniform along the corresponding wire sur- 
face. Cartesian {x, y, z) and cylindrical (p, ip, z) systems of coordinates are also 
shown with their origins lying on the waveguide axis. Curved arrow shows the 
azimuthal direction in which the polar angle ip grows. 

Po ~ 0, one recovers the solid cylinder of the radius d treated before [1] (note 
different length notations here and in Ref. [1]). Another asymptotics of the 
infinite annulus width, d = oo, transforms the system into the columnar defect 
of the radius po in the otherwise uniform medium that will be addressed below 
too. Cylindrical walls of the hollow waveguide support boundary condition, 
as described by Eq. (1) with the uniform along the length and circumference 
extrapolation lengths Ai„ for the inner surface and Aout for the outer one. We 
do not confine the values of Kin,out to be real concentrating on the properties 
of the structure at the complex de Gennes distances. For completeness, we also 
introduce the Aharonov-Bohm (AB) whisker [83-85] with its total magnetic 
flux $ coinciding with the cylinder axis. Similar to the solid cylinder case [1], 
our analysis will be based on linear Helmholtz equation (9) what, in the case of 
superconductors, together with the fact that the background magnetic field B is 
uniform means that we restrict our consideration to the range of the magnetic 
intensities and temperatures close to the transition between superconducting 
and normal states [31] even though the results obtained are covered by the full 
GL theory (see Ref. [1] for more discussion) and have much wider validity range. 
Geometry of the system dictates a natural choice of the cylindrical system 
of coordinates r = (p, cp, z) with its origin lying at the circles center and the z 
axis being parallel to the waveguide. We will seek the solutions of Eq. (9) with 



the vector potential written in the symmetrical gauge, 

2 
Then, one has 



X r 



TTp^ 



(17) 



A. (o, is, + 1^,0). (18) 

We will operate with the energy E through which the wave vector k is 
expressed, according to (4). Such a treatment describes a superconductor wire. 
A transition to the frequency uj of the acoustical or electromagnetic oscillations 
can be readily done with the help of (3). It is convenient at this point to choose 
dimensionless variables; namely, we will measure all lengths in units of the width 
of the ring d; accordingly, if not stated otherwise, all energies will be measured 
in units of the ground-state energy tt^It? /{2mp(P) of the infinite Dirichlet ID 
quantum well of width d; all momenta, in units l/d; magnetic fields, in units 
of h/{\q\(P)] magnetization, in units of h\q\/{2mp); 2D current density, in units 
of qh/ {nipd'^); current, in units of qh/{mpd); time, in units of 2mpd'^ /{ir'^h); 
magnetic flux, in units h/\q\. Discussion on the choice of the units and its 
relation to the description of the processes in different physical systems can be 
found in Ref. [1]. Then, Eq. (4) in the chosen units transforms to fc = tt\/E, 
and Eq. (9) for the order parameter '^{^]p,ip,z) becomes: 

1 a / a*\ 192^ 2i/i„ $\ a* /i„ $\\ 9^* 



-pd-p ['^) ^7W^-p U""" + ~p) ^" U""" + ~p) *+ a;^+"^* = ' 

(19) 
with E being a total energy of the particle. Factoring out the z-dependence 

*(<i>;p,^,z) = e^'=^^V('i>;P,^) (20) 

leads to the equation for the transverse function tp{'^; p, (f): 

p dp \ dp J p^ difi"^ P\2 P J dip \2 p J 

(21) 
Longitudinal wave vector k^ and the transverse energy E-^ are related as 



k^ = ny/E-E-^. (22) 

Rotational symmetry of the system allows one to separate out the angular and 
radial dependencies in the transverse function tp{^; p, (p): 

V'„™($;p,(^)-^e""'^i?„™($;p). (23) 

V 27r 



Here, m = 0, ±1, ±2, . . . and n — 0, 1, . . . are the azimuthal and the principal 
quantum numbers, respectively. In this way one arrives at the equation for the 
radial function _R„m($;p): 

d^Rnrn '^ dRnm ( "^ + ^ , ^ u \ 1} 2 



dp"^ p dp \ P 2 



Bp] Rnm+TT^Ei„„^Rnm^O- (24) 



This equation is supplemented by the boundary conditions for the function 
Rmni^', p)- Our choice of the vector potential in the form of (18) drops it out 
from (12) which becomes: 



dRn 



dp 



T 



1 



A.„ 



-R-n 



= 0. 



(25) 



Eqs. (24) and (25) constitute the problem of finding the eigenfunctions 
Rnm{^\ p) and eigencnergics i?$.„„ of the 2D circular annulus with its circum- 
ferences supporting the boundary conditions with, in general, complex Kin^out- 

Analytical solution to Eq. (24) is expressed via the Kummer confluent hy- 
pergeometric functions M{a,b,x) and U(a,b,x) [86-88]: 



R„ 
B 



»(*;P) =7$;nmexp ( -jpM i—p^ 



|m«|/2 



fu {pi,Aout)M ( a„^, |m$| + 1, yP^ ) - /m {pi,Kut)U ( a„ 
with 



|TO$| 



Pi = 


= Po + 1, 




?n$ : 


= m + $, 






TO$ + m$ + 1 


TT^ E^;nrn 



and functions Jm{PtA) and fir{p,A) written as 
|to$| B 



Im {p, A) = 
/c/(p,A) 



P 

|m$| 



M\ a. 



B 



U a 



TO$ 



m$ 



B 



B , 

1,— P 
' 2^ 

B , 

M,-p^ 






(27) 
(28) 

(29) 



TO$ 



h BpU' I Qm^, |m$| 



Here, prime denotes a derivative of the function with respect to the last ar- 
gument and the real coefficient 7$;nm is determined from the normalization 
condition which is either of the form from (7) for superconductors or 



po + l 



\R-i>;nrn{p)\'^ Pdp = 1 



(31) 



for the electromagnetic and acoustical waveguides. 

The form of the function Rnm{^', p) from Eq. (26) automatically satisfies 
condition (25) at the outer surface, p — po + l. Imposing boundary requirement 
at the inner edge, p — po, one arrives at the following transcendental equation 
for determination of the energies i?$.„^ : 



(26) 



B . 
1,— P' 

B , 



(30a) 
(30b) 



/m (po, -Ai„) fu (pi, Aout) - fu (po, -Ai„) /m (pi, Aotit) = 0. 



(32) 



In the limit of the vanishing background magnetic field, B ^ 0, Eqs. (26), 
(30) and (32) transform to (33), (34) and (35), respectively: 



Rnm{^; P) = l-i>,nm fv {Pl, Kut) J\ni^\ \^ \j E^.^^.^pj - /,/ {pi, Kut) ^Im^l [t^ \J E^.^^.^pj 

(33) 



1 



1 



Iy {P, A) = n^Ei^^„^Y{„^^^ [^^jEi,,mP) + X^"*! [^^Ein^p) (34b) 

fj{PO,~-^in)fY{Pl,Kut) - Iy {pa,-^in) f.i {pi, Kut) =0, (35) 

where J,y(x) and Y^{x) are Bessel functions of the first kind of the order v [86]. 
As expected, in the limit of the large radius, po -^ oo, Eq. (35) transforms into 
the corresponding dependence of the straight waveguide [11]: 

A,:„ \ r— / 1 



cos7ry^+ ( — -p= 7ry^Ai„ j sin7ry^=0. (36) 



Note that the azimuthal quantum number m disappeared at the transition from 
(35) to (36) since for the straight duct only the transverse index n defines the 
quantization. 

So far, no use has been made of the complexity of the de Gennes distances: 
Eqs. (26), (29) - (36) are valid for either real or complex parameters A. Mis- 
cellaneous cases of the real Robin parameter and different combinations of the 
background B and the AB fields for the same 2D geometry have been addressed 
before [15,24,48,52, 53,84,85,89-96] in the framework of the formalism of linear 
equations (2) or (9); in particular, Eq. (35) for the Dirichlet boundary conditions 
was presented in [84,85,89] with its uniform magnetic counterpart for zero AB 
flux given in [92,94] while Eq. (32) without the AB field was derived in Refs. 
[48,52]. 

From Eqs. (32) and (35) it straightforwardly follows that for the complex 
extrapolation lengths the energies, in general, are complex too. 









2 
with real -E$.^„ and r$;„m, and the following property holds: 



E,S>-nm (Am, Ao„t) — E^.,^^ (Ai„, Ao„t) (38) 

with the overline denoting a complex conjugate value. The same is true for the 
corresponding radial functions too. Physical meaning of the complex energy 
E from Eq. (37) was provided before [1]; in particular, based on the standard 



quantum mechanical theory of scattering [97,98], it was stated that its imagi- 
nary part is inversely proportional to the lifetime t of the corresponding quasi 

bound state: 

1 

T<S>;nm = -p • (39) 

Also, since Eq. (24) is invariant under the simultaneous transformations m — > 
mil, "Jj^'l'^Fl, the following relations hold: 

E<i,-nm — -E$±l:n,m=Fli (40) 

what is a general periodicity property of the AB systems [84, 85, 99, 100]. 

We repeat once again that even though (32) and (35) were derived for the 
linearized GL equation, it can be shown that for the thin superconducting rings 
they follow also from the complete nonlinear GL theory. For doing this, one 
needs to employ the method of minimizing of the free energy F of the super- 
conducting state with respect to the coefficient linking the order parameter of 
the nonlinear equation with its linear counterpart, Eq. (23) [52,53]. 

All other physical characteristics of the structure are found from the eigen- 
functions ^($;r) and eigenenergies E^-nm- For example, in the case of super- 
conductors a magnetization operator 

M..A_1V (41) 

is used to calculate the magnetic moment Mz of the 2D ring: 

M, = (V'„™($;/9,(^)|M|V'„™($;p,^)). (42) 

One immediately gets: 

D /-PO + l 



Mz = -lm + - J |i?„™($; p)r P''dpj . (43) 

An expression for the superconductor current density j^ in our dimensionless 
units is written as [30,31]: 

j, = Im [f (r) V* (r)] + A* (r) * (r) . (44) 

Similar formula (with, of course, A = 0) can be used for the Poynting vector 
in electrodynamics [101] or for the sound energy density flux in acoustics [102]. 
For the 2D case without the z dependence it transforms to 



. _ 1 

J^~2^ 



dR\ f m,i> 1 



^-(^^j^^Hp +r^"^'^ 



i|2 



--l^ 



(45) 



with the unit orthogonal vectors Gp and e^p along the radial and azimuthal 
directions, respectively. Since its divergence is proportional to F 

divj, = J|i?|2F, (46) 

10 



the positive (negative) imaginary part of the energy means that the correspond- 
ing spatial point serves as a sink (source) [1]. In addition, Eq. (45) shows that 
the total current Jp through the circle of the radius p is: 



Jp = p dipjsp = pim f R— 

In turn, for the 3D wire the total longitudinal supercurrent Jz is given as 



Jz 



po + 1 /•2ir 

pdp / dfjsz = 'Re{kz) cxp [-2lm{kz)z] . 

pa Jo 



(47) 



(48) 



Eq. (7) shows that the longitudinal dependence of the density n,, is of the same 
form: 

ns = \Rip)fcxp[-2lm{kz)z]. (49) 

Expressions for Re(fcz) and l-m{kz) were derived and analyzed before [1], so, 
here we simply rewrite them pointing only to the fact that the imaginary part 
of the wave vector defining, according to (48), the longitudinal alteration of the 
current vanishes together with T: 



Re(fc,) 
Im(A:^) 





(50a) 



(50b) 



E- E. 



(r) 



+ (r*;nm/2) + 



E-E. 



(r) 



In other words, the lossless longitudinal flow is achieved when the total trans- 
verse energy is real, as expected. As it will be shown below, for the multiply 
connected geometry it takes place when the total 2D radial current through the 
inner circle Jp^ is equal to its counterpart through the outer border Jp-^ . 

3 Results and discussion 

In this section the outcome of the calculations based on the theory developed 
in Chapter 2 is presented and their mathematical and physical interpretations 
are given. If in the results below the index or subscript $ is dropped, it means 
that its particular value of $ = is used; otherwise, it takes any arbitrary 
magnitudes. 

3.1 Columnar defect 

Before discussing the annular geometry with the finite inner and outer radii, 
it is instructive to consider the case of the cavity in the otherwise uniform 
material. In terms of Fig. 1 this corresponds to the infinite width d of the 



11 



ring, d ^ oo. The only remaining distance po obviously becomes the natural 
unit of length; accordingly, all other physical quantities containing distances 
will be correspondingly rescaled. This convention will be assumed throughout 
the whole present subsection for either real or complex Robin parameters. As 
there is only one surface, in the following the subscript 'in' at the extrapolation 
length A will be dropped. For our treatment, such a geometry presents an 
interest since it is an asymptotic case of the strong magnetic field that pushes 
the particle closer to the origin and, thus, the influence of the outer cylindrical 
surface becomes negligible. However, properties of the system with the hole 
in the film are important by themself with possible application, first of all, in 
superconductivity and they were addressed by the number of theoretical [96, 
103-109] and experimental [105,107,110,111] studies limited, however, to the 
real de Gennes distances with its, primarily, infinite value (Neumann boundary 
condition). Here, we will present the results that escaped an attention of the 
previous researchers. We will also make a comparative analysis with its inverse 
counterpart of the solid cylinder. Below, the geometry of the columnar defect 
will be called the exterior configuration, and the confined disk - the interior one 
[49,50]. 

Obviously, for the antidot with unrestricted at infinity motion of the particle, 
the bound states for the 2D geometry exist only when the background magnetic 
field is not zero, B 7^ 0. Assuming the same gauge for the vector potential as 
before and taking into account asymptotic properties of the confluent hyperge- 
ometric functions at infinity [86,87], one writes the following expression for the 
radial part of the total wave function: 

P ,^ , f B ,\fB ,V"-l/%^/^ m^ + |m^| + l tt^ i?i_ 

Rnm['^;P) =7'I>;nmexp \-JP ) ( "^^^ ) ^1 ^ ^ ^— J"l$ 

(51) 
Robin demand at the boundary p = 1 leads to the following transcendental 
equation for the determination of the energies E^.^^^ : 

. B 1\ /m$ + |m<j,| + l 7r2£^$.„„ b\ 

"^*""y-aJ^(^ 2 t^'I"^*i + 1'yJ 

+i?t/'(^^^^^^^^^-^%,l-*l + l,f)=0, (52) 

where, instead of the derivative of the function U, one can use the confiuent 
hypergeometric function itself with the different parameters according to [86] : 

U'{a, b, z) = -aU{a + 1,6+1, z). 

Eq. (52) for $ = was derived before [107] (see also Rcf. [108] for the Neumann 
limit). Normalization condition, Eq. (5), turns to 

\Rnm{<i>;p)\^pdp^l, (53) 
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and expression for the magnetization reads: 

M, = -(m+jJ |i?„„($;p)|'p3^pV (54) 

Four numbered equations above are valid for the real as well as complex A. 

Using properties of the function U{a,b, x) [86] , one immediately sees from Eq. 
(52) that for the small magnetic field, i? — > 0, the energy spectrum transforms 
into the Landau levels [98, 112] disturbed by the AB fiux [84,85,99, 100]: 

This is in a sharp contrast with the interior problem where the limit of Eq. 
(55) is achieved for the large intensities B [1]. In the latter case, the strong 
fields push the charged carrier closer to the center of the disk and, due to the 
small magnetic radius rs = B^^^^, the walls of the dot only slightly perturb 
the Landau- AB states, Eq. (55). In turn, for the columnar defect in the small 
fields, the cavity radius is much smaller than the diverging distance rs and, 
thus, the particle with energy from Eq. (55) almost docs not 'sec' the antidot 
and does not 'feel' the perturbation caused by it. 

3.1.1 Real de Gennes distance 

Another asymptotic limit of Eq. (52) can be analytically derived for the small 
nonzero magnetic fields and the extrapolation length tending to zero from the 
left, A — > —0. Then, one gets for the principal eigenvalues, n = 0, of each fixed 
m: 

^*;0,m = -^, A ^ -0, B ^ 0. (56) 

The same dependence is obtained from equation 

^^/\EiZj^U\ {-^^/l^ioj) + jl\rn,\ {^^J\Eioj) = 0, (57) 

describing the motion of the particle with negative energy E^.q^ inside the field- 
free disk with negative extrapolation length A [11]. Here, I,^{x) is a modified 
Bessel function of the order ly [86]. In the same way, for the annulus it can 
be shown from Eq. (35) that for the two (generally different) infinitely small 
negative extrapolation lengths Ain,out the two lowest energies -E$;o.m and -B*;!.™ 
for each m are: 



^*;(0,l),m - ~-5T2 ' Aj„ ^ -0, Ao„t ^ -0. (58) 



In fact, results (56) and (58) are a particular case of more general property of 
the vanishingly small negative Robin parameter intensively studied recently by 
mathematicians [113-117] (see also [2, 10]). Systems with negative extrapolation 
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Figure 2: Energy Eqq (left panels) and corresponding magnetization M^ (right 
panels) of the columnar defect as functions of magnetic field B and positive (up- 
per panels) and negative (lower panels) extrapolation length A. Note different 
scales for the upper and lower panels. To better emphasize the obtained de- 
pendencies, the A- and B- axes in panel (c) have been reversed as compared to 
panel (a). For the same reason, the field axis in panel (d) grows in the opposite 
direction of panel (b) . 



lengths raise interesting theoretical questions about fundamentals of quantum 
mechanics; for example, for an arbitrarily shaped domain with general perfectly 
reflecting walls a generalized Heisenberg uncertainty relation was derived ear- 
lier this year [118]. As it was mentioned in the Introduction, experimentally 
such systems have been fabricated with the help of superconductors [32, 35, 36]. 
According to Eqs. (15) and (16), the limit of the negligibly small negative de 
Gennes distance can be achieved by applying to such structures appropriately 
directed electric field £ when the total extrapolation length Ajot approaches zero 
from the left. 

Fig. 2 shows Eqq as a function of the positive [panel (a)] and negative [panel 
(b)] de Gennes distance A and magnetic intensity B. It is seen that the energy 
monotonically increases with the field for all extrapolation lengths. The asymp- 
totic limit of the function U{a, b, x) for the large positive x and negative a such 
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that a; < 26 - 4a [8( 
U{a, b,x) — 2exp 
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K (261 -sin 26*) + - (59) 



with K — &/2 — a and cos^ = x/(4«;), leads, for the Dirichlet case, A = 0, to 
the following expression of the energies i?$.„„j in the strong fields: 



Eib-r 



I _ 1 
Ia=o ^ 4,^2 



3\ 67r 



"+4/ B 



2/3' 



B^ + ^B, B 



OO. 



(60) 



As Fig. 2 demonstrates, the interplay between the magnetic field and the Robin 
electrostatic potential leads to the about the same quadratic dependence at the 
strong fields for all extrapolation lengths while the A-depcndence in the same 
limit is noticeable for the small positive Robin distances only when 1/A in Eq. 
(32) is comparable to B/2, and for the larger magnitudes of the de Gennes 
distance the energy E^.^^ is almost A- independent. 

The onset of the peculiarity from Eq. (56) is clearly seen in panel (b). It 
also shows that with increasing intensity B the negative energy grows too, turns 
to zero at some field and unrestrictedly increases in the positive direction with 
further growth of B. To clarify the interaction between the magnetic field and 
negative Robin parameter, it is instructive to consider the extrapolation length 
^*-m ^^ which the energy is zero, E^.^ ^ = 0. It directly follows from Eq. (52) 
that it is 



A 



(0) 
$;m 



I I ^ 

|W$| - — 



1 + (to$ + |m$| + 1) 



t/( "*+l;"*l+^ +l,|m.i,| + 2,f 



^( "^+l"^l+\ |m^| 



1, 



(61) 



Utilizing asymptotic properties of the confluent hypergeometric function U{a, 6, x)^, 
one gets in the limiting cases of the small and strong fields: 



1 



A, 



(0) 



ln(B/8)+7 T' ™* - U 

-{\m^\ + B/2), m^^Q 
-{B/2 + m^ + l), 



B<1 
B>1. 



(62) 



Here, 7 — lim„^oo (X]fe=i i ~ ^'^'^) ~ 0.57721 ... is Euler's constant [86]. 

Fig. 3 shows Am dependence on B for several quantum numbers m. The 
states with A lying above (below) the corresponding curves possess negative 
(positive) energies. Increasing field shrinks the range of extrapolation lengths 
with the energies lying below zero: for the smaller magnitudes of the de Gennes 



^Equations (13.5.21) and (13.5.22) in Rcf. [86] arc wrong. The correct expressions can be 
obtained from Eqs. (8.11) and (8.10) of Ref. [88] for the Whittaker functions Mi^^^(x) and 
Wk,^ij,{x) and their relations to the functions M{a, b, x) and U{a, b, x). 

^We point out other typos in the reference literature; namely, Eqs. (13.5.9) in [86] and 
(6.8.5) in [87] are wrong. For the correct form, one needs in Eq. (6.8.5) of Ref. [87] to convert 
the rightmost negative sign into its positive counterpart. 
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Exterior m=0 

Exterior m=1 

Exterior m=-^ 

Exterior m-2 

Exterior m=-2 

Interior m-0 



10 



B 



Figure 3: Critical extrapolation length Am of the columnar defect as a function 
of the magnetic field B where the solid line is for m = 0, dotted line - for m = 1, 
dashed line - for to = —1, dash-dotted curve - for m = 2, and dash-dot-dotted 
curve - for m — —2. For comparison, Aq of the quantum disk is also shown by 
the thick solid line. 
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distance the larger magnetic intensities are needed to pull out the energy into the 
positive area. For comparison, we also plot in Fig. 3 the critical extrapolation 
length of the dot calculated from [1] 
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Eq. (62) for the interior problem reads 

-52/16, 
|m$| - 



(63) 



1 



A 



(0) 



m<[, = I 
R^5/2), -*^0 I ^«1 (64) 

-(B/2 + |m<j,|), B>1. 



Comparison of Eqs. (62) and (64) shows that for the principal state, m,$ = 
0, much larger magnetic fields are needed for the disk geometry to reach a 
positive energy. This can be understood after consideration of the factors which 
affect energy behavior. For either case, tending to zero negative Robin length 
pushes energy downwards. For the quantum dot, an increasing magnetic field 
is the only parameter that forces the energy to grow. For the columnar defect 
it is additionally aided by the cylindrical electrostatic potential of the hole. 
Combined effective potential of the magnetic squeezing and electric confinement 
gets narrower for the larger B, and so, their mutual effort in withstanding 
against the opposite trend of the de Gennes distance requires smaller fields to 
see the energy i?$;o.m positive. 

Right panels of Fig. 2 show magnetizations Mz corresponding to the energy 
of the state depicted in the left parts. Its dependence on the field shows com- 
pletely different behavior as compared to the interior problem (see, e.g.. Fig. 9 
in Rcf. [1]). It is known [119] that the magnetic moment and the real energy E 
are related as 



M, = -TT- 



.dE 



(65) 



Accordingly, at the strong intensities, the magnitude of Mz for the Dirichlet 
case, A = 0, grows almost linearly with the field: 
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(66) 
as it follows from Eq. (60). Similar to the energy, the A dependence of magne- 
tization is noticeable in this regime for the small Robin parameters only. For 
the small negative de Gennes distance, A — > —0, the lowest state, n = 0, of 
each azimuthal quantum number m is localized near the interface with the ra- 
dial component of its wave function closer and closer resembling the (5-function; 
accordingly, the magnetization turns to 



Mz = - 



|l. A--0^ 



(67) 



oo. 
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3.1.2 Complex Robin parameter 

It was shown in Ref. [1] that nonzero real part of the complex extrapolation 
length A quenches the resonant features of the system parameters dependence 
on the imaginary component of A. As the same is true for the Robin antidot 
too, below we will present the results for the purely imaginary Robin parameter 
A = iA. 

ir] 

Fig. 4 depicts Eqq and Too as functions of A and B. The dependence of the 
real parts of the total energy on the field and on the complex Robin parameter 
is quite similar to that of the system with positive real extrapolation length [cf. 
panels (a) of Figs. 2 and 4]; namely, in either case, they almost quadratically 
depend on the magnetic intensity in the high-field regime while the magnitude 
of the Robin parameter |A| alters the energy at high fields only when it is of 
the order |A| ~ 2/B. In turn, the imaginary part exhibits a sharp extremum 
as a function of the purely imaginary de Gennes distance. The same resonant 
dependence was observed for the interior problem too [1]. The origins of these 
resonances are of the same nature: systems with almost purely Dirichlet case at 
|A| ^ 1 and that close to its Neumann counterpart, |A| ^ 1, are characterized 
by the same sign of T and its linear dependence on |A| and 1/|A|, respectively; 
accordingly, interaction of these two asymptotics in the intermediate regime, 
|A| ^ 1, leads to a pronounced maximum with its magnitude being n and m 
dependent. The drastic difference between the exterior and interior configura- 
tions lies in the T-B dependence: if, for the quantum disk, the field quenches 
the resonance, then for the hole, as panel (b) of Fig. 4 vividly demonstrates, 
the extremum increases with B. This inverse behavior is explained by the local- 
ization properties of the magnetic field: its growth, in the case of the quantum 
dot, pushes the charged carrier closer to the axis and, thus, the influence of the 
Robin potential on the properties of the system diminishes [1], while for the 
columnar defect it moves the particle closer to the circular impenetrable barrier 
and in this way increases scattering. Calculations show that the maximum of 
the resonance for the strong B almost quadratically depends on the field what 
is consistent with the result of the zero de Gennes distance, Eq. (60). Another 
difference between the two configurations is explained by the similar reasoning; 
namely, for the interior disk the resonance increases with the radial n and az- 
imuthal \m\ quantum numbers [1] and for the perforated film at the fixed field 
the opposite is true since the radial wave function for the larger n and |?7i| is 
spread out wider from the system origin and, thus, its deformation by the hole 
(disk) is smaller (larger). 

3.2 Annulus with complex extrapolation lengths 

Similar to subsection 3.1.2, our analysis of the hollow waveguide will concentrate 
on the purely imaginary Robin parameters: Ai„ o„( = i\in,out, since in this case 
the obtained energy and current dependencies on A are the most conspicuous 
ones. First, the results are discussed for _B = with the emphasis on the asymp- 
totic cases and the conditions under which the energy -E$.„^ is real. Influence 



18 




Figure 4: (a) Real Eqq and (b) negative double imaginary Too parts of the total 



^00 

?± 



transverse energy Eqq as functions of A and B for the columnar defect. Note 
reversed A- and _B-axes directions in panels (a) and (b). 



of the background uniform magnetic field B will be presented separately in the 
second subsection. 

3.2.1 Field- free case 

We start our analysis from the study of the asymptotic cases of the large and 
small |A|; namely, similar to the solid cylinder case [1], the limiting cases of the 
small, |Ai„^o«t| ^ 1, and large, |Ai„^o«t| 3> 1, magnitudes of the purely imaginary 
de Gennes distances can be derived from Eq. (35) using the addition theorem 
for the Bessel functions and their asymptotic properties [86]. Expressions for 
the imaginary part of the energy T^.„m are provided below for all possible 
permutations of the asymptotics: 
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is a Kronnecker symbol for, in general, real v and 
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^u'^ix) == J\u\{pox)Y\^\{pix)-Y\^\{pQx)J\^\{pix) (69a) 

F^'^ix) - J|,|(pox)l^;|(pix)- 1^,1 (pox) J|;|(pix) (69b) 

F^^'ix) = j[^^{pox)Y\,\{pix)~Y(^^{pox)J\,\{pix) (69c) 

J-,^^(x) = J{^^{pox)Y{,^{p,x)-Y{^^{pox)J{,^{p,x) (69d) 

and xJ/„ with i, j running over the indices D and A^ is nth solution of the 
equation zeroing the corresponding J^^^(x): 






0, i,j = DoTN 



(70) 



with the assumption that x^^ = 0. In other words, x^^^ are the eigenvalues of 
the 2D boundary problem with the Dirichlet (D) and/or Neumann (N) require- 
ments at the corresponding edge of the annulus with the first (second) character 
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Figure 5: Negative double imaginary part Too of the total transverse energy 
Eqq as a function of the radius po and equal inner and outer purely imaginary 
distances A ~ Ai„ = Xout- 

of the superscript denoting inner (outer) circle. Their dependence on the radius 
po for t^ = was analysed before [120]. In the asymptotic limit of the large 
radius, pa — > od, Eqs. (68) transform to 
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(71c) 



derivable, of course, also from Eq. (36). In Eq. (71b) the lengths A< and A> are 
defined according to |A<| = min(|Ai„| , |Ao„t|), |A>| = max(|Ai„| , |Ao„t|)- 

Fig. 5 demonstrates a transformation of the imaginary part of the energy 
with changing inner radius and equal outer and inner extrapolation lengths, 
A = Xin = Xout- It is seen that with increasing the radius po the magnitude of 
the resonance grows and, for po ^ 2 it saturates to its value of the straight waveg- 
uide. Larger scattering for the straight channel can be seen from the comparison 
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between Eq. (71a) and the corresponding expression r$;„m = 4j? ^, A/tt^ for 
the sohd cyhnder at the small |A| [1]. Here, jy„ is nth root of equation Ju{x) = 
[86]. In this regime, for m = n = $ = the imaginary part of the energy for 
the straight channel is about 3.41 times larger. Physically, this difference is 
explained by the centrifugal forces acting in the curved sample [121] and be- 
ing absent for the straight film where both interfaces contribute equally to the 
reactive scattering while for the annulus the unequal distribution of the wave 
function along the radius leads to the increased influence of one boundary and 
the decreased contribution from the other surface with their total mutual effort 
being smaller than for the unbent structure. 

Next, let us address the important question about the conditions under which 
the 2D energy E^.^^ is real. Apart from the theoretical interest related, for 
example, to the actively developing field of T^T-symmetric quantum mechanics 
[122-124] , this problem is of a paramount practical applications since real energy 
means a lossless longitudinal current down the waveguide. As it follows from 
Eq. (35), the energy E^.^^, in order to be real, has to satisfy at the same time 
the following equations: 



^^n^Z^ UJEi) XoutT^f,^ i^jEi) = 0, (72a) 



defining different branches of its Ai„-, Xout- ^nd po"dependence. Since for the 
fixed dc Gcnnes distances and inner radius the rea/ energy E;^.^^ has to obey two 
equations (72) simultaneously, this happens only for some definite correlations 
between Ai„, Xout and po- For example, it is immediately seen from Eqs. (68) 
that in the asymptotic cases the energy is real when the following conditions 
hold: 



Airi-' rn* [•'^nn^n) ~ Ao«t-/ rn$ [•'^m^n) i \Xin.out\ *C 1, ( ( oaj 

^'m^nj ^inJ^m^ \^m^n) ^ ~ ^^m^ {^m^nj/ ^out, | Ai„ | <C 1, | Ao„f | > 1, (73b) 

•^m« {^Tn^n) / ^in "^ {X„i^n) ^outJ^m^ \^m^n) ^ I Am | > 1, | Ao„t | <C 1, (73c) 

PO /Ai„ — -pi Aout , TO$ — n — U, J 






In turn, in the limit of the large radius, po — > oo, Eqs. (72) transform to 

(Ai„ + Aout) cos Try _B,| = 0, (74a) 

(A„Ao„t7r2^^ + l)sin7ry^ = 0. (74b) 

They also follow straightforwardly from the eigenvalue equation of the straight 
waveguide with Robin boundary conditions, Eq. (36). Then, one has following 
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23 
Figure 6: (a) One of the surfaces in (po, Am,Aout) space on which the energy 

Eqq is reaL (b) Real energy £',^5 corresponding to panel (a) as a function of po 

and Xin- 



situations 

E^ = 2\2 ' n^,n ^ 1,2,. . ■ for Ao«t = -Ai„, (75a) 

£;,-[ = (n + 1/2)2, n- 0,1,... f^^ A„Ao„t = ^^ ^ ^ . . (75b) 

TT^ [n + 1/2) 

Spectrum from Eq. (75a) has been known before [125-129] from T'T-symmetric 
[122-124] study of the quantum systems. This allows to understand the features 
of one of the sheets of the real energy surface in (po, Ai„, Xout) coordinates shown 
in Fig. 6(a). Namely, at the large radius the ground state has, for the small 
|Ai„|, as it follows from (75a), the real energy E^ = 1 under the condition 
Aout = — Ai„. As the inner Robin parameter passes the value of |Ai„| — I/tt, 
the ground energy becomes Eq — 1/ (7rAi„) , and the same condition for the 
inner and outer Robin distances holds. This steep energy descent is clearly 
seen in panel (b) of Fig. 6 exhibiting corresponding real energies Eq^. As the 
inner de Gennes distance goes through |Ai„| — 2/7r, the ground energy becomes 
Eq ~ 1/4, and the relation between Ai„ and Xout reads: Xout — ^4/ [ir'^Xin)- 
In other words, the described evolution is a transformation from the purely 
Dirichlet boundary condition Ai„ = Xout — with the energy spectrum being 
squares of integers to its Dirichlet-Neumann counterpart 1/Ai„ = Xout = with 
the spectrum of squares of half-integers through such intermediate values of 
Xin and Xout which guarantee a zero value of the imaginary part of the energy. 
Thus, the extra (compared to the pure Dirichlet case) state with the energy 
l/(7rA)2 in Eq. (75a) plays the role of the bridge linking two real-energy levels 
with the different types of the requirements at the walls. As Eqs. (72) and (73) 
demonstrate and Fig. 6 depicts, the similar dependence persists for the finite po 
too, however, the coefficients between the inner and outer purely imaginary de 
Gennes distances are now radius dependent. Of course, for the solid cylinder, 
Po = 0, the real energies £^$.„„ = {j\m^\n/'^)'^ are possible only for Xout = 0. 

To further elaborate on this issue, in Fig. 7 the complex energy Eqq [panel 
(a)] and the corresponding radial function -Roo(p) [panel (b)] are shown as func- 
tions of the inner de Gennes distance Xin for the fixed outer extrapolation length 
Xout and radius po- For the given fixed parameters, the energy turns real at the 
two imaginary Robin lengths: Ai„ = -0.9285438 and A,„ = -0.2033145. The 
form of the real parts of the corresponding radial functions in panel (b) mani- 
fests that first of them is close to the Neumann-Dirichlct edge of the spectrum 
while the main contribution to the second one is from the pure Dirichlet con- 
figuration, as it was discussed in the previous paragraph. To turn F into zero, 
each of the functions is deformed compared to the corresponding type of the 
boundary requirements and, in addition, contains an admixture of the imagi- 
nary component. However, even though the energies for these two parameters 
are real, the transversal radial current through the ring exists. This directly 
follows from the form of the functions in panel (b) of Fig. 7 and expression (47) 
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Figure 7: (a) Complex energy _E(^ as a function of the inner purely imaginary 

de Gennes distance \in for po = 0.5 and Aom* = 0.2 where the dotted line depicts 

{y\ 25 

the real part of the energy E)^^ while tne solid line is for the double negative 

imaginary part. Dashed horizontal line denotes zero energy, (b) Real (lower at 
p « po + 1) and imaginary (upper at p w Po + l) parts of the total radial functions 
Roo{p) for po = 0.5, Xout — 0.2 and several Ai„ where the thick solid lines are for 
Am = 0.2, dotted lines - for Am = 0, dashed lines - for Ai„ = —0.2033145, dash- 
dotted Hues - for Am — —0.35, dash-dot-dotted lines - for Am = —0.9285438, 
and thin solid lines are for Am = -1-2. Note that real parts vanish identically 



at n : 



1 



for the total radial current which for our parameters reads: 

Jp„,, ==±/^|i?(po4)l' (76) 

with plus (minus) sign corresponding to the inner (outer) surface. Note that in 
general case of the complex A, the expressions for these currents do not depend 
on the sign of the real parts of the extrapolation lengths: 

Jp,,=±Po,i|i?(Po,i)l'^^^;^^^^^- (77) 

\^^in.out I 

The crucial point here is the fact that the equality of these two fluxes forces the 
energy to be real: 

ri; -, =0. (78) 

•■I Pa — I PI 

In fact, Eq. (78) is a reflection of the current conservation law when the incoming 
flux through one of the surfaces is equal to the outgoing one and, thus, the 
longitudinal current is not affected by the processes occurring in the cross section 
of the waveguide. Expression for the divergence, Eq. (46), tells that every spatial 
point, including that at the annulus inner and outer circumferences, absorbs and 
emits the same amount of the flux and, thus, the total flow through it is zero. 
The difference with the systems with real Robin parameters is in the fact that 
in the latter situation no current passes at all through the points lying on the 
boundary. It immediately follows from equations (77) and (78) that for the real 
energy the following condition holds: 

[Im(A,„)Im(A„„0]lr=o<0- (79) 



Current density patterns corresponding to the wave functions from Fig. 7(b) 
are shown in Fig. 8. It is seen from panel (a) that for the same imaginary 
extrapolation lengths (or, more generally, for the de Gennes distances with 
the same signs of their imaginary parts, Xin^out > 0) the currents through the 
inner and outer interfaces flow in the opposite radial directions; accordingly, the 
annulus acts as a sink (source) for positive (negative) A with nonzero divergence 
of the current, and, as it follows from Eq. (46), no real energies can be obtained. 
The same is true for the configuration when one of the surfaces supports real 
Robin requirement. For example, even though the divergence at the Dirichlet 
inner border in panel (b) is zero due to the vanishing function on it, the outer 
reactively scattering wall does not allow to suppress the imaginary part of the 
transverse energy Eqq. In turn, for the real energies [panels (c) and (e)] the 
fluxes through the inner and outer circle head in the same direction and, since 
they exactly compensate each other, Jp^ = J^j^ , one gets a lossless longitudinal 
current with zero divergence of its transverse counterpart. Any deviation from 
these two points results in unequal currents through the confining surfaces when, 
depending on the ratio between imaginary de Gennes lengths with the opposite 
signs, the ring acts as a sink [panel (d)] or a source [panel (f)]. 
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Figure 8: Current densities of the state with n = m = $ = for po = 0.5, 
Aout = 0.2i and (a) A,„ = 0.2i, (b) A„, = 0, (c) A„ = -0.2033145^, (d) 
Ain = -0.35?, (e) A„ = -0.9285438z, and (f) Ai„ = -1.2i. For each of the 
figures, the currents are normahzed with respect to their largest value. The 
longer arrows denote larger currents. Corresponding energies £'^5 are written 
near each pattern. 
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3.2.2 Interplay of the magnetic field and imaginary de Gennes dis- 
tance 

For the quantum dot, the increasing magnetic field pushes the charged carrier 
to the centre of the disk, i.e., away from the Robin boundary with the corre- 
sponding decrease of its impact on the properties of the system; as a result, 
the resonance on the F — Im(A) characteristics gets smaller with growing B 
and, ultimately, in the limit of the very strong intensities, one recovers the 
states with independent of A real energies in the form of the Landau levels, Eq. 
(55) [1]. For the ring geometry with nonzero po the formation of the Landau 
states is prevented by the inner impenetrable wall which, in the case of the 
growing fields, plays the more and more dominant role since the distribution 
of the parameter \['(r) concentrates around it. Accordingly, one can expect 
the enhancement of the reactive scattering for the nonzero value of the imagi- 
nary part of its Robin parameter. To exemplify this, we present in Fig. 9 the 
energies Ej^q of the annulus with radius po = 0.1 for the two configurations 
of the extrapolation length A: the left panels show real and imaginary parts 
of the total transverse energy for the opposite signs of the Robin parameters 
on the inner and outer surfaces, A — Xout — —Km while for the right panels 
the Dirichlet boundary condition was adopted on the inner wall: A = Xout, 

(r) 

Ain — 0. At the zero field, the real parts £"$.„„ of the energy show with 
the increasing magnitude of the de Gennes distance a transformation from the 

Dirichlet case at A = with i?$.„„ = i^u,^ \n/''^) to the either pure Neu- 
mann [panel (a)], i?^.„„ = (^i^^in/'"') 7 or Dirichlet-Neumann [panel (c)], 

^k-nm — \^Wn \nl''^] ' geometry at the large Robin parameter, \^ 1. Imagi- 
nary parts at -B = exhibit discussed above variation of F$;„m with its passage 
through zero for the opposite signs of A [panel (b)] or the resonant dependence 
caused by the outer imaginary extrapolation length [panel (d)]. Growing mag- 
netic field restricts the motion to the area around the inner surface with the 
corresponding change of the energies i?$.„„; for example, for the inner negative 
A one sees in panels (a) and (c) of Fig. 9 that, for the strong fields, the real part 
^i-nm varies with the de Gennes distance only for A < 2/ B, while the F depen- 
dence transforms into the resonance profile with almost quadratic dependence 
on _B - the behaviors familiar from the study of the columnar defect, section 
3.1.2. For the Dirichlet inner wall in the increasing field, the particle moves 
away from the outer scattering Robin interface with the corresponding decrease 
of its infiuence on the system; accordingly, the F dependence on A quenches and, 
at the large B, tends to zero for all external extrapolation lengths [panel (d)]. 
At the same time, the real part grows with the field with its dependence on A 
fading too, and in the limit of the strong fields such that Bp^ ^ 1, one gets: 
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Figure 9: Real [panels (a) and (c)] and double negative imaginary [panels (b) and 
(d)] parts of the total transverse energy Eqq as functions of the purely imaginary 
Robin parameter A and field B for the annulus with inner radius po = 0.1. For 
the left panels A = Xout = —Kn while for the right panels A = Xout with Am = 0. 
Note reversed A and B axes in panel (b) as compared to other panels. 



(recall that we use different length units here and in subsection 3.1). In our 
model of the electric field influence on superconductors, Eqs. (15) and (16), it 
means that the combined application of the external electric £ and magnetic B 
fields allows to vary the scattering in a broad range switching from the lossless 
current to any desired degree of scattering by the simple change of the gate 
voltage and/or magnetic intensity. Comparing inner Dirichlet case of the ring 
with the energy E^.^^ dependence on Im(A) and B for the disk (Fig. 2 in Ref. 
[1]), one sees their qualitative and, for the imaginary part, almost quantitative, 

(r) 

similarity with the difference being in the E^.^^ dependence on the field for the 
large B: linear, Eq. (55), for the disk and quadratic, Eq. (80), for the annulus. 
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4 Concluding remarks 

Ubiquitousness in Nature of the physical systems which are described by the 
wave equation (2) supplemented by the Robin boundary condition, Eq. (1), 
with complex extrapolation length A dictates a necessity of their correct theo- 
retical description. In the present research, the previous results for the simply 
connected domain of the 2D circular disk and 3D solid cylinder [1] were supple- 
mented by the calculations of the simplest doubly connected geometry of the 2D 
ring and infinitely long 3D cylinder of the annular cross section with, in general, 
different complex Robin lengths on the two confining circumferences. Com- 
parative analysis between the two configurations revealed new features for the 
multiply connected structures; for example, it was predicted that under some 
correlation between the outer and inner Robin distances with the opposite signs 
of their imaginary parts, the eigenvalue of Eq. (2) turns real - despite of the 
complexness of the boundary condition (1). Physically, this is explained by the 
same amount of the fiux entering the system through one interface and leaking 
out from it through the other one. This result, interesting by itself for, e.g., 
the 7'T-symmetric quantum mechanics [122-124], is of a large practical signif- 
icance too since the real transverse energy £^$.„^ restores lossless longitudinal 
current down the waveguide. New phenomena are predicted to exist for the 
system response to the external magnetic field B too; namely, if, the growing 
B for the disk decreases the imaginary part of the complex transverse energy 
turning, for the large intensities, the levels into the Landau states, then, for the 
2D ring, F fades with the increasing field only if the inner Robin parameter is 
real, otherwise, due to the reactive scattering at the inner wall, the imaginary 
part of the energy grows quadratically with B. 

We considered the hollow structures invariant under the rotation around the 
z axis. In the process of the sample growth intentional or unintentional devia- 
tions from the azimuthal symmetry can occur. Such Neumann structures have 
been studied both theoretically [130-132] and experimentally [133,134]. The 
same Dirichlet geometry was used for the study of the distorted semiconductor 
rings [135, 136]. Combined infiuence of the excentricity, real de Gennes distance 
and uniform magnetic field B was also calculated [55]. Based on the results 
presented above, the conjecture can be made that if the displacement psh of the 
inner disk from the centre of the outer one is larger than its radius po, then 
for the high enough fields when the ratio rs/ipsh ^ Po) is smaller than unity, 
the imaginary part of the energy will asymptotically tend to zero with its real 
part approaching the Landau levels, Eq. (55). In the opposite case, psh < Po, 
the T — B characteristics at the strong fields will depend on the imaginary part 
of the inner extrapolation length Ai„; e.g., it will decrease to zero for i? — > oo 
at Im(Ai„) — and its magnitude will almost quadratically increase for any 
nonzero Im(Ai„) while the real part of the energy E^^' in either case will ex- 
hibit the same ~ B^ dependence. The waveguides with more complicated cross 
sections (with, for example, two or more leaking holes inside the larger circle 
with complex Robin parameter) can be treated in a way being a combination of 
the methods developed here and in the cited references, and, depending on the 
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particular geometry, they should exhibit one of the dependencies just described. 
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